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Abstract 

Stars closely approaching the massive black hole in the center of the Galaxy provide a unique opportunity to 
probe post-Newtonian physics in a yet unexplored regime of celestial mechanics. Recent advances in infrared 
stellar spectroscopy allow the precise measurement of stellar Doppler shift curves and thereby the detection of 
(3 2 post-Newtonian effects (gravitational redshift in the black hole's potential and the transverse Doppler shift). 
We formulate a detection procedure in terms of a simplified post-Newtonian parametrization. We then use 
simulations to show that these effects can be decisively detected with existing instruments after about a decade 
of observations. We find that neglecting these effects can lead to statistically significant systematic errors in 
the derived black hole mass and distance. 

Subject headings: black hole physics — gravitation — relativity — techniques: radial velocities — Galaxy: 
center — infrared: stars 



1. INTRODUCTION 

General Relativity (GR) is the least tested theory of the four 
fundamental forces of nature. The m ~ 3-4 x 10 6 Mq mas- 
sive black hole (MBH) in the Galactic center (GC) and the 
stars around it (Eisenhauer et al. 2005; Ghez et al. 2005) pro- 
vide a unique laboratory for detecting post-Newtonian (PN) 
effects and for probing GR (for a review, see Alexander 2005). 
Long-term astrometric and spectroscopic monitoring of stars 
orbiting the MBH allow the derivation of their orbital ele- 
ments. The stars are observed at periapse distances as small 
as r p 10 3 r s (r s = 2Gm/c 2 ) with velocities as high as 
f3 p ~ few x 0.01 (J3 = v/c). The relativistic parameter at pe- 
riapse, T = r s /r p > 10 -3 is a few times higher than that on 
the surface of a white dwarf. This is an unexplored regime of 
celestial mechanics. For comparison, the high precision con- 
firmation of GR predictions in the Hulse-Taylor binary pulsar 
(PSR 1913+16) was for masses of ~1.4M S with (3 P ~ 0.003 
at r p ~ 2xl0 5 r s (Taylor & Weisberg 1989). Note that although 
some accretion processes occur at T > O(0. 1), the complex 
physics of accretion and its many uncertainties severely limit 
the usefulness of accretion emission as a probe of GR. 

At present all the observed orbits can be adequately mod- 
eled in terms of Newtonian physics. With improved resolu- 
tion, higher precision and a longer baseline, deviations from 
Keplerian orbits may be detectable. In this Letter we examine 
the future detectability of leading order ( ft 2 ) PN effects with 
existing instrumental capabilities. We show that such effects 
can be detected in the observed stellar Doppler shift curves 
(note that the Doppler shift is no longer equivalent to the ra- 
dial velocity when PN effects are included). 

The measurement of stellar proper motions and accelera- 
tions in the inner arcsecond (0.04 pc) of the GC (Eckart & 
Genzel 1996; Ghez et al. 1998) preceded by more than a 
decade the measurement of radial velocities for these stars 
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(Ghez et al. 2003; Eisenhauer et al. 2003a). Thus, most the- 
oretical studies of the detectability of PN orbital effects fo- 
cused on those that affect the proper motions, such as the 
periapse shift (e.g. Rubilar & Eckart 2001; Weinberg et al. 
2005). However, with the advent of adaptive-optics-assisted 
IR imaging spectroscopy (e.g. VLT/SINFONI; Eisenhauer et 
al. 2003b), it is now the observed Doppler shifts rather than 
the proper motions that provide the tightest constraints on de- 
viations from Newtonian orbits. The current quality of stel- 
lar spectroscopy in the GC allows the determination of the 
Doppler shift, z, to within 5z ^25 kms _1 /c, or SX/X^ 10~ 4 
(Eisenhauer et al. 2005). 

The observed Doppler shift curve z(t) can be expanded in 
terms of the 3D source velocity (3{t'), where t 1 and t are the 
light emission and arrival times, respectively, 



z = AA/A = B + fli/3 + B 2 f3 2 + C?(/3 3 



(1) 



As shown in §2, two effects contribute to the 2nd-order term 
(Kopeikin & Ozernoy 1999), B 2 = B 2 ,t + B 2 , g : the special 
relativistic transverse Doppler effect, with B 2yt — 1/2, and the 
GR gravitational redshift, with B 2 . g = 1/2. It then follows 
that the 2nd-order PN effects are detectable with existing in- 
struments, since (B 2 , t + B 2 ^)/3j-10" 3 >(5A/A-10" 4 . 

This Letter is organized as follows. In §2 we discuss the 
properties of the (3 2 terms and the implications for their de- 
tectability. In §3 we demonstrate the detectability of these 
effects by simulated observations. We discuss our results in 
§4. 

2. f3 2 EFFECTS IN THE DOPPLER SHIFT 

When T 1 the orbit is Keplerian to a good approximation 
and the orbital energy equation (the vis viva equation) is 



f ~r s /r-r s /2a, 



(2) 



where a is the orbital semi-major axis. For an eccentric orbit 
r p «a and so near periapse T ~ 1 . Two effects contribute to 
the (3 2 term of the Doppler shift expansion (Eq. 1). One is the 
GR gravitational redshift, 



z g = (l- r s /ry 1/2 - 1 ~ r s /2 



(3) 



which also contributes a constant to the 0th order term of the 
z-expansion, z g = r s /4a+(l/2)/3 2 = B 0tg + B 2 ^/3 2 (Eq. 2). 
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FIG. 1 . — Bottom: the full relativistic radial velocity curve of S2 near peri- 
apse. Top: the contribution of the post-Newtonian ffl effects of gravitational 
redshift and transverse Doppler shift to the total. 

The second is the full relativistic Doppler shift of a moving 
source seen by an observer at rest, 

z D = (l+/3cosi?) I Vl-/3 2 - 1 ~/3costf + /3 2 /2, (4) 

where § is the angle between the velocity vector and the line 
of sight. The Newtonian Doppler shift, z r = (3 cos d — Bi(3, is 
the 1st- order term in the z-expansion. The special-relativistic 
transverse Doppler effect, z t = /3 2 /2 = B 2 .tP 2 , contributes to 
the 2nd-order term in the z-expansion. 

The light travel time from the star to the observer changes 
with the orbital phase in an inclined orbit. We include 
this effect, known as the R0mer delay, by relating the ob- 
served Doppler shift z(t) to the actual velocity at emission 
/3{t—[Ro — r(t') coailj(t')]/c}, is the angle between r and 
the line of sight and i? is the distance to the MBH. It is cru- 
cial to do so because the phase dependence of the light travel 
time also contributes to the 2nd order term (e.g. Alexander 
2005). 

The gravitational redshift and the transverse Doppler shift 
affect only the measured Doppler shifts, but not the astromet- 
ric proper motions. These are treated here as Newtonian, since 
astrometric PN effects are below the current astrometric pre- 
cision. Figure 1 shows the 2nd-order PN effects that are pre- 
dicted for the star S2, based on the orbital elements derived for 
it from a Newtonian fit to the data (Eisenhauer et al. 2005). As 
expected, the PN effects near periapse are a few percents of 
the total measured Doppler shift, since for z t + z g <C 1, 

(z t + z g )/z~ 0$/0 p ~ fa. (5) 

Brumberg et al. (1975) studied the detectability of the grav- 
itational redshift and transverse Doppler shift effects in binary 
pulsar timing data. They show that radial velocity data alone 
are not enough to detect these effects, since to leading order 
they can be absorbed in a purely Newtonian solution by very 
small modifications of the orbital elements K (radial velocity 
amplitude), lu (argument of the periapse) and Vb (the line- 
of-sight center-of-mass velocity). To break the degeneracy 
they suggest using the GR periapse shift (w ^ 0). While this 
may be feasible for a binary pulsar, where the orbital period 
is very short and the timing is extremely accurate, a different 
approach is required for the orbits around the Galactic MBH. 
In this context there are two possibilities to break the degen- 
eracy. One is to use the astrometric data, which can constrain 



both lu and K. Another is to use the fact that the stars all or- 
bit a common center-of-mass and share the same Vo- Thus, 
a combined solution of the astrometric and radial-velocity or- 
bits for all the stars can break the degeneracy and reveal the 
2nd-order PN effects. 

Best-fit Newtonian orbital elements and their errors are usu- 
ally obtained by a x 2 fit of Newtonian orbits to the data. We 
extend this procedure by introducing a new parameter, X 2 , 
which quantifies the strength of the 2nd-order PN effects, so 
that the GR gravitational redshift term is parametrized as 

z g (X 2 ) = (l-X 2 r s / r y 1/2 -1, (6) 
and the full relativistic Doppler term as 

z D (X 2 ) = (l + 0cos0) /^l-X 2 ^ -1. (7) 

For Newtonian physics, X 2 = 0, whereas for PN GR, X 2 = 1. 
Values other than zero or one are also possible in principle and 
would indicate deviations from the predictions of both theo- 
ries. The X 2 parametrization is chosen here for its simplicity, 
but is neither general nor unique. For example, it arbitrar- 
ily presupposes that deviations from GR (expressed by z g ) 
are related in a specific way to deviations from special rela- 
tivity (expressed by zd)- If future data quality or candidate 
alternative theories warrant it, the fit procedure can be easily 
generalized to include multiple PN parameters. 

We looked for evidence of PN effects in the existing orbital 
data by applying the PN fit procedure to the orbits of the 6 
stars SI, S2, S8, S12, S13 and S14 (Eisenhauer et al. 2005). 
We obtain a best-fit value of X 2 = —4.9 ±2.9, which appears 
to exclude PN GR at a 0.02 confidence level and Newtonian 
physics at a 0.05 confidence level. The number of currently 
available Doppler shift measurements is very small, and the 
astrometric data is not homogeneous. It is conceivable that 
these results are strongly influenced by systematic effects. Fu- 
ture observations will settle this question. 

3. SIMULATIONS 

We performed a suite of simulations to investigate the de- 
tectability of 2nd-order PN effects in the Doppler shift data 
with existing instrumental capabilities and with realistic astro- 
metric and spectroscopic precision. Since the relative strength 
of the PN effects scales with (3 p (Eq. 5), we focus on stars that 
pass through deep (high velocity) periapse during the moni- 
toring period. This has actually happened in the period 2000- 
2003 (Eisenhauer et al. 2005; Fig. 2) for the two stars S2 
(with period P = 15.2 yr, eccentricity e = 0.88, r p = 1500r s , 
P P = 0.02) and S14 (with P = 38 yr, e = 0.94, r p = 1400r s , 
P p = 0.03). Unfortunately, at that time precise stellar spec- 
troscopy of the inner GC was just becoming available, and 
only a few Doppler measurements of S2 were obtained (Ghez 
et al. 2003; Eisenhauer et al. 2003a). 

To make the simulations realistic, we used the 6 orbits 
that were actually monitored and solved by Eisenhauer et al. 
(2005). We generated synthetic orbits based on the published 
Newtonian orbital elements. The simulated datasets consisted 
of one observation per month (2 astrometric coordinates and 
one Doppler shift per star), for 6 consecutive months per year 
(the GC is observable only part of the year). The starting point 
was the year 2000, so the first 3 years already cover the peri- 
apse passage of S2 and S 14 (Fig. 2). To each measurement we 
added uncorrected random Gaussian errors, with a standard 
deviation of 25kms _1 for the radial velocities and 1.5 mas 
per coordinate for the astrometry, which are representative of 



PROBING POST-NEWTONIAN GRAVITY 



3 




0.4 0.3 0.2 0.1 -0.1 2000 2002 2004 2006 2008 2010 

Aa [arcsec] Time [years] 



FIG. 2. — Simulated orbits of the stars S2 and SI 4. The "observed" astro- 
metric data (left) and Doppler shift data (right) are marked by points. 



the measurement errors for S2. The MBH's astrometric po- 
sition and proper motion were assumed to be known, since 
frequent accretion flares can reveal the MBH's location (Gen- 
zel et al. 2003b; Ghez et al. 2004). These assumptions are 
somewhat optimistic, since S2 is better measured than most 
stars, and since accurate localization of the MBH by the flares 
is still challenging. However, judging by the continuing im- 
provement in the data acquisition, it is not unreasonable to 
expect such data quality in the near future. 

We generated progressively longer datasets, starting from 
3 years and up to 18 years. We repeated the entire proce- 
dure twice, both with and without the PN effects. We then 
applied the PN orbital fit procedure to the simulated datasets, 
and recorded the best-fit values of X 2 , as well as those of m 
and R . 

The upper panel of Figure 3 shows the derived values of 
X 2 for both Newtonian and PN simulated orbits of the 6 stars 
(of which only S2 and S14 go through deep periapse). The 
errors on X 2 clearly improve with longer monitoring. The 
two hypotheses can be decisively distinguished by the data 
at a 4cr confidence level after 10 years of monitoring. The 
lower panel translates the errors on X 2 into expected detec- 
tion significance for 1 (S2 only), 2 (S2 and S14), and all 6 
stars. These confidence levels reflect the assumed measure- 
ment errors. The effect of a change in the assumed errors can 
be easily assessed by noting that a global scaling of all errors 
simply scales the confidence interval by the same factor. 

4. DISCUSSION 

We demonstrated that leading order PN effects in the 
Doppler shifts of stars orbiting the Galactic MBH can be de- 
cisively detected in about a decade of astrometric and spec- 
troscopic monitoring with existing instruments. We added to 
the usual set of Keplerian orbital elements a simple PN pa- 
rameter, X 2 , which quantifies the strength of the PN effects. 
This PN parameter can be easily generalized to a set of PN 
parameters to describe any extension of Newtonian physics. 

It is important to optimize the observing strategy to improve 
the chances of detecting PN effects. Generally, the observa- 
tional program should focus on deep periapse passages and 
include as many stellar orbits as possible to better constrain 
the parameters and break the Newtonian/PN degeneracy. Our 
simulations were all based on the currently observed orbits 
and assumed optimal observing conditions. However, future 
observations will include additional orbits and be carried un- 
der varying observing conditions, which were not considered 
here. A comprehensive study of a range of orbital configura- 
tions, measurement errors and observational resources alloca- 
tion strategies remains to be done. 

Radial velocity information allows the unambiguous 
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FIG. 3. — Top: Ict limits on the best-fit X2 parameter that quantifies the 
strength of the leading order PN effects in the data (X2 = for Newtonian 
physics and X2 = 1 for PN physics). Bottom: The mean significance of 
detection of the PN effects (rejection of Newtonian hypothesis) in the simu- 
lated dataset with PN physics, as function of monitoring time for 1, 2 and 6 
stars, presented in terms of the rejection probability (left axis) and as standard 
deviations (right axis). 



derivation of both the MBH mass, m, and distance to the 
GC, Ro, without prior assumptions (Salim & Gould 1999). 
If ignored, the PN effects can introduce systematic errors in 
the best-fit values of these parameters. For example, the in- 
clusion of the gravitational redshift, transverse Doppler and 
R0mer effects (not taken into account in previous orbital so- 
lutions) in the orbital fitting of existing data changes the 
best-fit values from m = (3.47 ± 0.25) x 10 6 M Q to m = 
(3.65 ± 0.28) x 10 6 M Q and from R Q = 7A7 ± 0.24 kpc to 
R = 7.59 ± 0.25 kpc. These changes are not yet statisti- 
cally significant. However, our simulations indicate that once 
Doppler shift measurements are available, the discrepancies 
can become highly significant. For example, the neglect of 
the PN effects in the orbital solutions of the simulated data of 
6 stars monitored for a decade results in 8er discrepancies in 
the derived values of m and Rq . 

Another closely related post-Newtonian effect is phase- 
dependent stellar flux variability. This is due to the com- 
bined effects of relativistic beaming and the Doppler shift. To 
leading order, beaming causes the bolometric flux to scale as 
1-4/3 cos i? (Rybicki & Lightman 1979, Eq. 4.97b). However, 
since the infrared bands lie in the Rayleigh-Jeans part of typi- 
cal S-star spectra, the flux in the emitted spectral range that is 
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observed after being Doppler-shifted scales as 1 + 3/3 cost?. 
Therefore the total variability scales as 1 — (3 cos d, which 
translates to an expected variability near periapse of order 
/3 P ~ 0(0.01). 

We note that the effects of any dark distributed mass around 
the MBH on the Doppler shift is likely to be small, but per- 
haps not completely negligible. The change in the Doppler 
shift at periapse, Sz^m, due to the retrograde periapse shift, 
8oj dra , that is induced by the gravitational potential perturba- 
tion of the dark mass (see Alexander 2005) is of the order 
fcd m ~ 0(|<5wdm| P P ). This perturbation can be neglected as 
long as l^dml *C $ P ~ 0.01. The periapse shift due to the 
hypothesized particle dark matter cusp is completely negligi- 
ble (e.g. Gnedin & Primack 2004). The periapse shift in the 



orbit of S2 due to the observed luminous stellar mass (Genzel 
et al. 2003a) is only |<5wdm| ~3x 10~ 4 . The possible shift due 
to an extreme case of mass segregation of stellar black holes 
(Alexander & Livio 2004) may be as high as \Swd 10" d . 

To summarize, we have shown that recent advances in in- 
frared stellar spectroscopy allow the decisive detection of /3 2 
post-Newtonian effects in the Doppler shift curves with exist- 
ing instruments after about a decade of observations. 
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